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We find that the spin-charge motion in a strongly confined two-dimensional hole gas (2DHG) 
supports a propagating mode of cubic dispersion apart from the diffusive mode due to momentum 
scattering. Propagating modes seem to be a generic property of systems with spin-orbit coupling. 
Through a rigorous Keldysh approach, we obtain the transport equations for the 2DHG, we analyze 
the behavior of the hole spin relaxation time, the diffusion coefficients, and the spin-charge coupled 
motion. 



The correct description of the motion of charge and 
spin packets in materials such as semiconductors requires 
appropriate transport equations. While transport equa- 
tions for the Fermi liquid usually describe separate diffu- 
sion processes for spin and charge densities [l| , spin-orbit 
(SO) coupling alters their form dramatically. For exam- 
ple, the presence of coupling between spin and charge, a 
direct result of SO coupling, leads to the possibility of 
controlling the spin-orientation by electric fields 2] . An- 
other example is the recently proposed spin-charge mode 
propagating with a velocity equal to the Rashba SO 
coupling. This mode has no counterpart in usual trans- 
port and is intimately related to a non-equilibrium ver- 
sion of the spin-current in the ground-state of the system 



In the case of the two-dimensional electron gas (2DEG) 
with Rashba spin-orbit interaction, diffusive transport 
equations are known 0, 0. A similar analysis for the 
two-dimensional hole gas (2DHG), where the spin-orbit 
coupling is much stronger, has so far been lacking. In this 
letter we obtain the transport equations for very thin p- 
doped quantum wells in which only the heavy-hole band 
is populated. We start with a general Keldysh equation 
for spin-1/2 systems with generic SO coupling, and then 
particularize to the case of the heavy hole 2DHG. The 
same approach can be applied other types of SO cou- 
plings, such as the 3D Dresselhauss with little or no mod- 
ification. The equations obtained here form the basis for 
studies of transport and spin-orbit effects in 2DHG sys- 
tems. We find that, unlike in the case of 2DEG, the spin 
orbit effects do not produce first order derivative terms in 
the transport equations, but appear only in second and 
even third order derivative terms, which now have to be 
taken into consideration. As such, the 2DHG does not 
exhibit the uniform spin polarization under an applied 
electric field, or the spin-galvanic effect that its 2DEG 
counterpart does. The diffusion coefficient becomes a 
tensor with components for spin and charge which de- 
pend on the spin-orbit coupling. Its spin component 
turns out to usually be smaller than the charge compo- 
nent, which might provide an alternate explanation for 



the recent experimental observations pj attributed to the 
spin-coulomb drag. The Dyakonov-Perel relaxation time 
for the heavy holes' spin is calculated as well. We dis- 
cuss the relation of the previously calculated spin Hall 
conductivity in the ballistic regime || E| and this for- 
malism. In addition, the transport equations support 
a spin-charge propagating mode, analogous to the one 
predicted in Ref. 3]. However, its dispersion relation is 
changed because of the different form of the SO inter- 
action, and is cubic rather than linear in momentum. 
The effects of Coulomb repulsion on the spectrum of the 
propagating mode are also accounted for in the present 
treatment. 

The most general single electron Hamiltonian with 
spin-orbit coupling for spin 1/2 systems is: 
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with \i(k) = — Xi(— k) due to time reversal invar iance. 
The energy eigenvalues are E±(k) = k 2 /2m ± |A(fc)| 
where |A| = \fX[\i. The retarded and advanced Green's 
functions of the clean 2DHG have the form : 
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where 
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k 2 1 (2m) and e is measured 
from the Fermi energy cp — k 2 F j2m. The density matrix 
g = g(k, r, t) is a 2 x 2 matrix in spin indices, whose inte- 
gral over the momentum space yields the reduced density 
matrix p = p(r, t) : 
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-g(k, r, t) = n(r) + Si(r)vi (3) 



where n(r, t), Si(r, t) are the charge and spin densities, d 
is the space-dimension, and v is the <i-dimensional density 
of states. The charge-spin density g satisfies a Keldysh- 
type equation which takes electron scattering into ac- 
count through an isotropic momentum relaxation time 
r: 
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This quantum Boltzmann equation can be solved iter- 
atively in the case of small spatial gradients Vg of the 
charge-spin density @|. Fourier transforming with re- 
spect to time and defining a linear functional K(g) = 

^{G R p — pG A ) — \ |f^, Jt~}> one obtains an equation 

for g(k,r,cj): 

. (2A 2 - n 2 )K(g) + 2X i X j a i K(g)a j 



continuity equation expanded up to the third order in 
spatial derivatives has the form: 
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flXlcr K(a)] where z, j, k = {x,y}.We first concentrate on the coeffi- 
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with = uj + i/r. Expanding in gradients g = g° 
SgW + 5g( 2 ) + . . ., a perturbative chain of equations is 
obtained up to the desired accuracy. Integrating over 
the momentum and taking the DC limit cor << 1 one 
obtains the transport equations for p. In the case of the 
Rashba coupling A^ = ottijkj this was done in Ref. p. 

For 2D systems with an isotropic Fermi surface, which 
includes the Rashba model, the 2DHG model, and the 
2D Dresselhauss model, in the limit of epr <C 1, e a = 
|A(fcp)| <C ep the continuity equation takes the form: 
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where 5g^ contains first order spatial derivatives, and 
subsequent terms contain higher order derivatives. Here 
m = 1, 2; m/3. The first term on the right hand side 
has zero components in the charge-sector and represents 
spin-relaxation with characteristic time: 
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where £ = e a r. Equation ® and the relaxation time for- 
mula J7J) are valid for any spin-orbit coupling in a 2D 
system with an isotropic Fermi surface. The general ex- 
pressions for higher order derivative terms tend to be 
exceedingly long. 

We now turn to a particular case: the two-dimensional 
hole gas. The valence band of type III — V bulk semi- 
conductors is a spin 3/2 band, separated into heavy and 
light-hole bands which touch at the T point. In a quan- 
tum well a gap opens at the T point so that the light 
hole band is moved below the heavy hole band. A suf- 
ficiently thin well ensures a gap large enough to place 
the light-hole band below the Fermi level. Then only the 
heavy-hole Hamiltonian should be considered. Asymme- 
try of the quantum well introduces a spin-orbit term of 
the form A^ with: 



X x = ak y (3k x —k y ), X y = ak x (3k y — k x ) (8) 

and X z = 0[l0|. This model has an isotropic Fermi sur- 
face with |A| = ak 3 . The limit of validity of our calcula- 
tions is h/r, ak 3 F <C kp/2rrihh where rrihh is the heavy- 
hole mass. For a compact notation, we define the charge- 
spin four- vector = (n, 5), \i = 0, 1, 2, 3. The generic 



cients Cf v and . In the 2DHG case we find Cf v 







due to the symmetries of the Hamiltonian. This is in con- 
trast with the linear Rashba model of the 2DEG, where 
nonzero linear terms were calculated in Ref. 3] . The dif- 
fusion coefficient tensor has only three distinct ele- 
ments due to the x y symmetry of the 2DHG Hamil- 
tonian. As we have already mentioned the general forms 
for the transport coefficients are very cumbersome so we 
assume small SO coupling (e a <C ep) and expand up to 
linear order in ^ in all coefficients (D's and i^'s) except 
for the charge diffusion coefficient whose first correction 
due to SO coupling is at order (f^) 2 . We find 
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where Do = -§-, and with symmetries — > 

D yy, D l\ -> D fy Another notable symmetry is = 
—D%% so the off-diagonal terms in the sum D^didjN u 
cancel. The expressions above are valid in the limit when 
the spin-orbit band splitting and momentum relaxation 
energy are assumed to be smaller than the Fermi energy 
i.e. e a , - « ep. However, the value of £ is unrestricted 
and we find for some common physical parameters that 
< x < 1. The transport equation above with these dif- 
fusion coefficients represents one of the main results of 
this paper. 

Here we briefly show how an interesting new phe- 
nomenon, the spin Hall effect, emerges from the above 
equations. Under an applied electric field along the x- 
axis we expect out-of-plane spins to be transported in 
the y direct ion [HI H^. We take into account the effects 
of an electric field per the substitution Vn — > Vn — evE 
0. From the transport equations above we see that E x 
generates: 



J' = {Df x ev)E- 
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which is in the preferred formal and where the sym- 
metry D£y = —Dyx gives us the antisymmetric tensor 
structure J\ = crsH^ijxE x . In the low impurity limit i.e. 



FIG. 1: (color online) (a) Charge diffusion coefficient (b) Spin 
relaxation time plotted vs SO coupling a(c) S z diffusion co- 
efficient (c) S y diffusion coefficient. All 3d figures are plotted 
vs. SO coupling strength a and hole- doping 5. 



as r — > oo we find a value of 
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This gives us: J^J = . When we take into account, 

like |8| that we are using spin- 1/2 matrices to represent 
a spin-3/2 system, we account for an additional factor of 
3 which yields: 



9e 



(13) 



in agreement with |8|, |9| . This value for the conductiv- 
ity matches the clean-limit and supports the previously 
known fact that vertex corrections vanish for this model. 







Nontrivial dynamics, in the absence of applied fields, 
appears in the third order derivative terms of the trans- 
port equations. The sum H^ k didjdkN v only depends 
on coefficients H££ x , H^y y1 and symmetrized combina- 
tions: m v 
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Hyyx + ^xyy • The number of distinct coefficients is fur- 
ther reduced to just four by the exact symmetries of the 
Hamiltonian: 
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FIG. 2: (color online) Third order coefficients plotted vs. SO 
coupling strength a and hole- doping 5. 



For the third-order coefficients we find: 
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where Ho = and up to this order we have the addi- 
tional symmetries: H xxx = ^xxx an d ^xxx = ~^xxx- 

It was found in Ref.j^ that in 2DEG systems the 
SO interaction can produce a mode of coupled spin and 
charge oscillations (fractionalized spin packets mode, or 
spin-galvanic mode) which can be weakly damped even 
in the diffusive regime. It was later found that the 
Coulomb interactions of the charge component do not 
suppress this mode Here we study the analogue 

of such mode in the 2DHG system using the ansatz 
5n(x,t) = 5nexp[iqx — iut], S y {x,t) = S y exp[iqx — iut], 
S x , S z = to find a solution of Eqs. (0 in the absence of 
external fields. Following the procedure of Ref. 13] to ac- 
count for Coulomb interactions, we obtain the spectrum 
of the mode 
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where k = ql s , l s = \^Dor s is the spin-diffusion length, 
R = 1+2ttkb/k and Kb = ve 2 l s characterize the strength 
of the Coulomb interaction, and 7 = H xxx r s /l 3 . As in 
the 2DEG case, the shape of the spectrum is controlled by 
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FIG. 3: Real (solid line) and imaginary (dashed lines) parts 
of the mode frequencies u)±(k) in an infinite 2DHG structure 
with Coulomb repulsion. Two windows of Re(co) / exist. 
The grey areas are outside the limits of the diffusive approx- 
imation. 



two dimensionless parameters 7 and Kb- The estimate for 
7 ~ € a T p /kFls can be compared with the 2DEG estimate 
J2DEG ~ £^ DEG T s /kFls- One observes that the larger 
SO splitting e a associated with 2DHG systems can be 
compensated by the small factor r p /r s . The Coulomb 
interaction parameter is large, Kb ~ 10 3 for l s = lOOnra. 

A sketch of the spectrum is shown in Fig. 03 As in 
the case of 2DEG, the modes are purely dissipative at 
k = 0, representing the charge diffusion and spin relax- 
ation. Spin and charge motions get coupled at n > and 
two intervals with Re(uo) 7^ emerge. As in the case of 
2DEG, the position of these intervals is very sensitive to 
the strength of the Coulomb interaction Kb- However, in 
addition to that, the Re{uj)/ Im(uo) ratio is also affected 
by the Coulomb interaction in the 2DHG case. Our cal- 
culation shows that the mode is always highly damped in 
the first interval. The mode can be weakly damped with 
Re(oo) / Im{uj) ^> 1 in the second window n > ft*, however 
ft* ~ (kb/'J 2 ) 1 ^ 3 is usually too close to the boundary of 
the validity of the diffusive approximation which requires 
n < ^t s /t p . We conclude that observing the fractional- 
ized spin packets propagating mode in the 2DHG requires 
more fine tuning of parameters than in the 2DEG case. 
However, the presence of the weakly damped solution at 
the boundary between the diffusive and ballistic regimes 
suggests that a search for a ballistic counterpart of the 
mode has to be performed. 

In conclusion, in the present letter we derived the dif- 
fusive transport equations for the 2DHG system. Our 



calculation shows that the components of the diffusion 
coefficient tensor depend on the magnitude of the SO cou- 
pling and that the spin diffusion is generally smaller than 
the charge diffusion, thus offering an alternate explana- 
tion to spin-coulomb drag for this already experimen- 
tally observed 7] phenomenon. We find that, in contrast 
to the 2DEG case, the spin-charge coupling enters into 
the transport equations starting with the second order 
derivative terms. Our diffusive equations reproduce the 
value of the spin Hall coefficient, previously calculated for 
this system in the ballistic regime H E| , thus supporting 
the conclusion about the absence of vertex corrections in 
this 2DHG system. The fractionalized-spin packets prop- 
agating mode in the 2DHG is also studied. It is found 
that the Coulomb interaction suppresses this mode more 
than in the case of electron gas. The spectrum of the 
mode suggests that its counterpart may exist for 2DHG 
systems in the ballistic regime. 
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